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Abstract 

Low-frequency quadrupole vibrational modes in deformed 36,38,40]Yjg close to the 
neutron drip line are studied by means of the quasiparticle-random-phase approxi- 
mation based on the coordinate-space Hartree-Fock-Bogoliubov formalism. Strongly 
collective K'^ = 0"*" and 2+ excitation modes carrying 10 — 20 Weisskopf units in 
the intrinsic isoscalar quadrupole transition strengths are obtained at about 3 MeV. 
There are two reasons for the enhancement of the transition strengths. First, the 
quasiparticle wave functions generating these modes possess spatially very extended 
structure. The asymptotic selection rules characterizing the (3 and 7 vibrations in 
stable deformed nuclei are thus strongly violated. Second, the dynamic pairing ef- 
fects act strongly to enhance the collectivity of these modes. It is suggested that the 
lowest K'^ = 0^ collective mode is a particularly sensitive indicator of the nature 
of pairing correlations in deformed nuclei close to the neutron drip line. 

Key words: Hartree-Fock-Bogoliubov method, Quasiparticle-RPA, Collective 
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1 Introduction 

The physics of drip-line nuclei is one of the current frontiers in nuclear struc- 
ture physics [1,2,3]. The number of unstable nuclei experimentally acces- 
sible will remarkably increase when the next generation of radioactive ion 
beam facilities start running. We shall be able to study the properties not 
only of the ground states but also of low-lying excited states of drip-line 
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nuclei in the medium-mass region. Collective excitation in neutron-rich nu- 
clei is one of the most interesting issues in this field. Because properties 
of low-frequency collective vibrational modes are quite sensitive to surface 
effects and details of shell structure, we expect that new kinds of collec- 
tive excitations emerge under such new situations of nuclear structure. In 
order to quest for collective modes of excitation unique to Tinstablc nuclei 
associated with new features such as neutron skins, many attempts have 
been made using the self-consistent RPA based on the Skyrme-Hartree-Fock 
(SHF) method [4,5,6] and the Quasiparticle-RPA (QRPA) including the pair- 
ing correlations [7,8,9,10,11,12]. A number of similar approaches using differ- 
ent mean fields have also been carried out [14,15,16,17,18,19,20,21,22]. (Sec 
Refs. [12,19,23] for extensive lists of references concerning the self-consistent 
RPA and mean-field calculations.) Most of these calculations, however, are 
restricted to spherical nuclei. 

Quite recently, low-frequency RPA modes in deformed nuclei close to the neu- 
tron drip line have been investigated by several authors. The time-dependent 
Hartree-Fock method formulated in the three-dimensional coordinate space 
with a complex absorbing boundary condition was applied to low-frequency 
isovector dipole modes [24]. Possible appearance of low- frequency octupole 
vibrations built on superdeformed states in neutron drip-line nuclei was dis- 
cussed in Ref. [25] on the basis of the SHF plus mixed representation RPA [26,27,28] 
calculations. In Ref. [29], we investigated properties of octupole excitations 
built on superdeformed states in neutron-rich sulfur isotopes by means of the 
RPA based on the deformed Woods-Saxon (WS) potential in the coordinate- 
space mesh-representation. We found that low-lying states created by excita- 
tion of a single neutron from a loosely bound low-il state to a high-f2 resonance 
state (Q being the z-component of the angular momentum) acquire extremely 
strong octupole transition strengths due to very extended spatial structure of 
particle-hole wave functions. All of these calculations, however, did not take 
into account the pairing correlation. In Refs. [30,31], low-lying Gamow- Teller 
/9-decay strengths were investigated by means of the proton-neutron RPA us- 
ing the SHF + BCS approximation. Gamma vibrations in ^^Mg were studied 
using the QRPA with the BGS approximation [32] on the basis of the response 
function formalism. It should be noted that these calculations relied on the 
BGS approximation, which is inappropriate, because of the unphysical nucleon 
gas problem [34], for describing continuum coupling effects in drip line nuclei. 

The nature of pairing correlations in neutron drip-line nuclei is one of the 
most important subjects in the physics of unstable nuclei. One of the unique 
features of drip-line nuclei is that the pairing correlation takes place not only 
among bound levels but also including continuum states. To describe this 
unique character of pairing, the coordinate-space Hartree-Fock-Bogoliubov 
(HFB) formalism is suitable [33,34] and has been widely used for the study 
of single-particle motion and shell structure near the continuum [35,36,37,38]. 
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Due to the pairing and continuum effects, spatial structure of quasiparticle 
wave functions near tlie cliemical potential changes significantly, which affects 
the properties of low- frequency excitation modes [39]. In order to study the ef- 
fects of pairing on the low-frequency excitation modes in deformed nuclei near 
the neutron drip-line, we have extended the previous work to self-consistently 
include pairing correlations, and constructed a new computer code that carries 
out the deformed QRPA calculation on the basis of the coordinate-space HFB 
formalism. 

The aim of this paper is to carry out the deformed QRPA calculation for 
neutron drip-line nuclei and investigate the low-frequency quadrupole vibra- 
tional modes with K'^ = O"*" and 2"*" in 36,38,40]yjg ^iQgg ^j-^g neutron drip line. 
According to the Skyrme-HFB calculations [40,41] and Gogny-HFB calcula- 
tion [42], these isotopes are well deformed. The shell- model calculation [43] 
also suggests that the ground state of ^°Mg is dominated by the neutron 
two-particle-two-hole components, which is consistent with the breaking of 
the N = 28 shell closure discussed in [44]. We investigate properties of low- 
frequency modes of excitation in these Mg isotopes simultaneously taking into 
account the deformed mean-field effects, the pairing correlations, and excita- 
tions into the continuum. 

This paper is organized as follows: In the next section, the framework of the 
mean-field and QRPA calculations is briefly described. In Section 3, results 
of the RPA calculation for low-frequency quadrupole vibrations with — 
0+ and 2+ in 36,38,40]y[g ^^^.^ presented and discussed focusing our attention 
to the microscopic mechanism of emergence of collective modes in deformed 
superfluid nuclei close to the neutron drip line. Concluding remarks are given 
in §4. 

A preliminary version of this work was previously reported in Ref . [45] . 



2 Method 



2.1 Mean-field calculation 



In order to discuss simultaneously effects of nuclear deformation and pairing 
correlations including the continuum, we solve the HFB equation [33,34,46] 



h''{ra) - h^{ra) 



\ 



[ }f{ra) -(h-(ra)-X-)) \^ljra) 



(1) 
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directly in the cylindrical coordinate space assuming axial and reflection sym- 
metry. In comparison to the conventional method of using a deformed har- 
monic oscillator basis, this method is more effective in the treatment of spa- 
tially extended wave functions, hke loosely bound states, resonant states and 
continuum states. As is well known, when the quasiparticle energy E is greater 
than the absolute magnitude |A| of the chemical potential, the upper compo- 
nent ifiira) obeys the scattering- wave boundary condition, while the lower 
component <^2(^"cr) is always exponentially decaying at infinity. 

For the mean-field Hamiltonian /i, we employ the deformed Woods-Saxon po- 
tential with the parameters used in [29] , except the isovector potential strength 
for which a slightly smaller value, 30 MeV in stead of 33 MeV, is adopted in 
order to describe "^^Mg as a drip-line nucleus in accordance with the Skyrme- 
HFB [40,41] and Gogny-HFB calculations [42]. The pairing field is treated 
self-consistently by using the density-dependent contact interaction [47,48], 



1 - P 



T) 

V ^"0 



S{r-r'), (2) 



with Vq = —450 MeV-fm^ and Qq = 0.16 fm^'^, where g^^{r) denotes the 
isoscalar density and P^- is the spin exchange operator. For the parameter rj, 
which represents density dependence, we use rj — 1.0 (surface type). Sensitivity 
of calculated results to the parameter r) will be examined in subsection 3.4. 
The pairing Hamiltonian is then given by 



It/ \ ^0 

h (r) = — 
^ ^ 2 




f{r). (3) 



^'0 

The normal and abnormal (pairing) densities are given by 

^^(P,^) = E E \v>UP^z,a)\', (4) 

a a=±l/2 

Q^{P,Z) = -Y1 E ^l,a{P^Z^'^)^2,aiP,Z,<^) (5) 
a a=±l/2 

and the mean-square radii of protons and neutrons are calculated as 

2x ^ J pdpdzr^g^{p,z) 
^''^^ Jpdpdzgr{p,z) ^ 

where r = {p, z) , r = Vp^ + and t—tt or i/; g^{p, z) and Q^i^p, z) being the 
proton and neutron densities. The average gaps are defined by [34,37] 

(A,) = - j dr^^(r)/i^(r)/ j drg^{r). (7) 



We construct the discretized Hamiltonian matrix by use of the finite differ- 
ence method for derivatives and then diagonalize the matrix to obtain the 
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quasiparticle wave functions on the two-dimensional lattice consisting of the 
cylindrical coordinates p and z. The kinetic energy term and the spin-orbit 
potential are evaluated using the 9-point formula. Because the time-reversal 
symmetry and the reflection symmetry with respect to the x — y plane are as- 
sumed, we have only to solve for positive Q and positive z. We use the lattice 
mesh size Ap = A2; = 0.8 fm and the box boundary condition at pmax = 10.0 
fm and ^max = 12.8 fm. The quasiparticle energy is cut off at 50 McV and the 
quasiparticle states up to fi'^ = 13/2''^ are included. This model space is larger 
than that used in Ref. [45]. We impose the condition on the convergence of the 
pairing energy as \ {Epl^^ — E^'l'^) / Epl^^\ < 10~^, where i denotes the iteration 
number and the pairing energy is defined by [37] 

We use the same deformation parameter P2 — 0.28 in the Woods-Saxon poten- 
tial for both neutrons and protons. This parameter is chosen to approximately 
reproduce the Q— moments calculated in Ref. [40]. We checked that properties 
of the QRPA modes do not change significantly when the deformation param- 
eter is varied around (^2 ~ 0.3. 



(8) 



2.2 Quasiparticle- RPA calculation 



Using the quasiparticle basis obtained in the previous subsection, we solve the 
QRPA equation in the standard matrix formulation [49] 



E 

7<5 




(9) 



This method is convenient to analyze microscopic structures of the QRPA 
eigenmodes in comparison with other RPA formalisms based on the Greens 
function method. 



The residual interactions in the particle-particle channel appearing in the 
QRPA matrices A and B are self-consistently treated using the density-dependent 
contact interaction (2). On the other hand, for residual interactions in the 
particle-hole channel, we use the Skyrme-type interaction [50] 



Vph[r,r 



5(r-r'), (10) 



with to = -1100 MeV-fm^, = 16000 MeV-fm^xo = 0.5, and 2:3 = 1.0. 
Because the deformed Wood-Saxon potential is used for the mean-field, we 
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renormalizc the rcsidiial interaction in the particle-hole channel by multiply- 
ing a factor /p/j to get the spurious K'^ = 1+ mode (representing the rota- 
tional mode) at zero energy {vph —>■ fph ■ Vph) ■ We cut the model space for the 
QRPA calculation by + < 30MeV, which is smaller than that for the 
HFB calculation. Accordingly, we need another self-consistency factor fpp for 
the particle-particle channel. We determine this factor such that the spurious 
= O"*" mode associated with the number fluctuation appears at zero energy 
{vpp —>■ fpp ■ Vpp). In the present calculation, the dimension of QRPA matrix is 
about 3700 for the K"" = 0+ modes in ^°Mg. 

In terms of the nucleon annihilation and creation operators in the coordinate 
representation, iplra) and ip'^lra), the quadrupole operator is represented as 
Q2K = J2a I drr'^Y2K{r)ip^ {ra)ip{ra) . The intrinsic matrix elements {X\Q2k\^) 
of the quadrupole operator between the excited state |A) and the ground state 
|0) are given by 

(Aig^i^io) = E Q^2KUf^, + g^p) = E m^kU. (n) 

a/3 a(3 

where 

<52K,l/3 = 27r5x,n„+n^ jdpdzQ^2K,ap (p^ ^) > (12) 

with 

<92K,l^(P, Z) = pWl,a{p, Z, i)v'2,/3(P, Z, f) - V?l,a(p, Z, T)V'2,/3(P, Z, j) 

- <^1,/3(P, Z, i)(^2,a(p, Z, t) + (pi,f}{p, Z, T)(^2,a(p, Z, i)}Q2K{p, z). (13) 

Here Q2k{p, z) = Q2K(r)e-''"^ = r-'Y^Kie, ^)e-''"^. 
We calculate the transition strength functions 

s'^ico) = j:\{\\QfKm'5ihw - hw,) (14) 

A 

for isoscalar quadrupole operators Q^j^ = Q2K + Q2KJ use notations 
B{Q'^2) — |(A|(55i<:|0)P for transition strengths and Mj. — {X\QIk\^) for tran- 
sition matrix elements (r = tt, i^, IS). Note that these quantities are defined 
in the intrinsic coordinate frame associated with the deformed mean field, so 
that appropriate Clebsh-Gordan coefficients should be multiplied to obtain 
transition probabilities in the laboratory frame [51]. For instance, a factor 1/5 
should be multiplied for obtaining the transition strength B{E2; 2f 0^) 
from the 2i state to the 0^ state, while the factor is unity for obtaining the 
transition strength B{E2; Ogg — * 2^) from the ground state to the 2^ state 
built on the excited K'^ = 0+ state. Here, 2'1 denotes the 2+ member of the 
ground-state rotational band, while 0^ and 2^ indicate the rotational band 
members associated with the K'^ — 0'^ intrinsic excitations. 
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Fig. 1. Single-particle energies in the deformed WS potential for neutrons in ^*^Mg, 
plotted as functions of the quadrupole deformation parameter P2- Solid and dotted 
lines denote positive- and negative-parity levels, respectively. Single-particle levels 
are labeled with the asymptotic quantum numbers [TVnsAjfi. 



3 Results and Discussion 



3. 1 Some features of calculated results 



The single-particle shell structure around the Fermi surface for neutrons in 
36,38,40]yjg exhibits an interesting feature. Figure 1 shows the single-particle 
energy diagram for the WS potential as functions of deformation parame- 
ter (32- As increases, a level crossing between the up-sloping [303] 7/2 level 
and the down-sloping [310] 1/2 level takes place, and a deformed shell gap is 
formed at A^" = 28 around (32 = 0.3. This deformed closed shell approximately 
corresponds to the (/7/2)~^(P3/2)^ configuration in the spherical shell model 
representation. The highest occupied level in this deformed closed shell is sit- 
uated very near to the continuum threshold, so that there is no bound level 
above it. However, neutron particle-hole excitations may take place into res- 
onance levels hke [303]7/2, [301]l/2 [312]3/2 lying just above the continuum 
threshold. In fact, as we shall discuss below, these resonance levels partic- 
ipate in the pairing correlations and play an important role in generating 
low-frequency collective modes of excitation in 36,38,40]y[g^ Thus, '^°Mg and its 
neighboring isotopes provide an interesting situation to investigate collective 
modes unique in unstable nuclei near the neutron drip line. The resonance 



Table 1 

Ground state properties of ^^'^^'^'^Mg obtained by the deformed WS-HFB calculation 
with /?2 = 0.28. Chemical potentials, average pairing gaps, and root-mean-square 
radii for protons and neutrons are listed. 







A. 




A. 


A. 




nucleus 


(MeV) 


(MeV) 


(fm) 


(MeV) 


(MeV) 


(fm) 


36Mg 


-20.0 


0.0 


3.06 


-2.09 


1.81 


3.74 


38]V[g 


-23.0 


0.0 


3.08 


-1.15 


1.98 


3.86 


40Mg 


-25.1 


0.0 


3.10 


-0.41 


2.14 


3.99 



character of these levels just above the continuum threshold is confirmed by 
means of the eigenphase-sum method (see Appendix) . 

Results of the deformed WS plus HFB calculation for the ground state proper- 
ties of 36,38,40 jyjg ^^.^ lig^g^j Table 1. Calculated values of the average pairing 
gap for neutrons are rather close to the value estimated in terms of the con- 
ventional systematics Agyst = 12 /^/A ~ 1.9MeV. On the other hand, the 
average pairing gaps for protons vanish. As shown in this table, the neutron 
root-mean-square radius increases as approaching the neutron drip line, while 
the proton root- mean-square radius remains almost constant. This means that 
the neutron skin structure emerges in these nuclei; the difference between the 
neutron and proton radii in ^'^Mg is about 0.9 fm. 

Results of the QRPA calculation for quadrupole transition strengths are dis- 
played in Fig. 2. We see prominent peaks at about 3 MeV for both the K'^ — 0+ 
and 2"*" excitations. Their strengths are much larger than the single-particle 
strengths indicating collective character of these excitations. The strength of 
the lowest K"" = 2+ excitation gradually increases as approaching the neutron 
drip hne, while the lowest K'^ — 0+ excitations in ^^Mg and ^°Mg seem to be 
split into two peaks in the case of ^^Mg. In the following, we make an extensive 
analysis on microscopic structure of these low-frequency collective excitations. 

3.2 K"" ^ 0+ modes 

We first discuss the K'^ = 0+ excitation modes in "^°Mg. The QRPA transi- 
tion strengths are compared with unperturbed two-quasiparticle strengths in 
Fig. 3 . A prominent peak is seen at about 3.2 MeV in the QRPA strength 
distribution; it possesses an enhanced strength of about 22 Weisskopf unit (1 
W.u. ^8.1 fm*^ for "^'^Mg). From the QRPA amplitudes listed in Table 2, it is 
clear that this collective mode is generated by coherent superposition of neu- 
tron excitations of both particle-hole and particle-particle types. In Fig. 3, the 
QRPA strengths are also compared with the strengths without the dynamical 
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Fig. 2. Isoscalar quadrupole transition strengths B{Q^^2) for the K = 0'^ excita- 
tions (upper panel) and the if = 2+ excitations (lower panel) built on the prolately 
deformed ground states of ^^'^^'^OMg. The arrows beside the abscissa axes indicate 
the neutron threshold energies, -Eth = 4.06 MeV (one-quasiparticle (Iqp) contin- 
uum; |A| +miiiEa), 4.12 MeV (two quasiparticle (2qp) continuum; 2|A|) for ^^Mg, 
2.31 MeV (2qp continuum) for ^^Mg and 0.82 MeV (2qp continuum) for ^^Mg. 
The QRPA calculations are made by using the surface-type pairing interaction and 
P2 = 0.28 for both protons and neutrons. 



Table 2 

QRPA amplitudes of the K"" = 0+ mode at 3.2 MeV in ^OMg. This mode has 
B{E2) = 3.4 e^fm^, 5(Q''2)=136 fm^, and B{Q^^2) = 182 fm^. The single-particle 
levels are labeled with the asymptotic quantum numbers [A^risAJfi of the dominant 
components of the wave functions. Only components with — Iffa/?]^ > 0.01 are 

listed. 





a 


/3 


(MeV) 




^20,a/3 

(fm^) 


2U,ap 

(fm2) 


(a) 


i/[310]l/2 


i/[310]l/2 


3.54 


0.438 


6.36 


4.27 


(b) 


z/[301]l/2 


i/[310]l/2 


3.93 


0.067 


-2.57 


0.925 


(c) 


i/[312]3/2 


i/[312]3/2 


3.99 


0.280 


-2.03 


1.08 


(d) 


z/[301]l/2 


i/[301]l/2 


4.32 


0.027 


0.992 


-0.176 


(e) 


z/[303]7/2 


i/[303]7/2 


5.76 


0.077 


-3.39 


0.966 


(f) 


z/[321]3/2 


i/[321]3/2 


7.15 


0.011 


3.23 


0.396 
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Fig. 3. Left: Isoscalar quadrupole transition strengths B{Q^^2) for the K'^ = 0+ 
excitations in ^^Mg. Results of the QRPA calculation with and without including 

the dynamical pairing effects are plotted in the upper and middle panels, respec- 
tively, while unperturbed two-quasiparticle strengths are shown in the lower panel. 
Notice that different scale is used for the unperturbed strengths. The arrow beside 
the abscissa axis indicates the neutron threshold energy 2|A| = 0.82 MeV. Right: 
Two-quasiparticle excitations generating the lowest = 0"*" mode at 3.2 MeV. The 
single-particle levels for the deformed WS potential are labeled with the asymptotic 
quantum numbers [NnsA]Q,. The chemical potential A is indicated by the dashed 
line. 



pairing effects, i.e., the result of QRPA calculation ignoring the residual pair- 
ing interactions. One immediately notice that the transition strength to the 
lowest excited state is drastically reduced when the dynamical pairing effects 
are ignored. 

Let us discuss the reason why the lowest = 0^ mode acquires eminently 
large transition strength. There are two points to understand this mechanism: 
1) existence of unperturbed two-quasiparticle configurations possessing large 
transition strengths, and 2) effect of residual interactions producing coherence 
among various two-quasiparticle configurations. 

To examine the first point, we plot in Fig. 4 spatial distributions of the 
quadrupole transition amplitudes for major two-quasiparticle configurations 
generating the lowest K'^ = 0+ mode. We see that they are notably extended 

beyond the half-density radius. This is a situation similar to that encountered 
in Ref. [29], where a neutron excitation from a loosely bound state to a reso- 
nance state brings about very large transition strength. We also note that the 
transition strength associated with the z/[301]l/2 ® z/[310]l/2 configuration is 
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'r"-l"""" ' i ' 

10 02468 10 024 
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Fig. 4. Spatial distribution functions Q^2^ap^P^ ^"^^ some two-quasiparticle excita- 
tions generating the lowest K'^ = 0+ mode in ^°Mg. The contour lines are plotted at 
intervals of 0.002. The solid and dashed lines represent positive and negative quan- 
tities, respectively. The thick solid line indicates the neutron half-density radius; 
^^(0)/2 ~ 0.045fm-3. 

much enhanced although it should be hindered if the selection rule AA/" = 2 
for the asymptotic quantum numbers is applied. This selection rule is broken 
for matrix elements associated with loosely bound states, because their radial 
wave functions are spatially extended and quite different from those of the the 
harmonic oscillator potential. 



Concerning the second point, wc have found that the dynamical pairing plays 
an especially important role. This point is easily seen by comparing the QRPA 
calculations with and without the dynamical pairing effects, which are shown 
in Fig. 3. It is apparent that the prominent lowest peak disappears when the 
dynamical paring effects are ignored. We can say that the coherent super- 
positions among the particle-hole, particle-particle and hole-hole excitations 
are indispensable for the emergence of this mode. The importance of the cou- 
pling between the (particle-hole type) j3 vibration and the (particle-particle 
and hole-hole type) pairing vibration has been well known in stable deformed 
nuclei [51]. A new feature of the K'^ — 0"^ mode in neutron drip-line nuclei 
under discussion is that this coupling takes place among two-quasiparticle 
configurations that are loosely bound or resonances, so that their transition 
strengths are strikingly enhanced. In addition, as seen in Fig. 4, their spatial 
structures (peak positions and distribution) are rather similar with each other. 
This is a favorable situation to generates coherence among them [39] . The im- 
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Fig. 5. Left: Isoscalar quadrupole transition strengths B{Q 2) for tlie 
K'^ = 0+ excitations in ^^Mg are plotted in the upper panel, while unperturbed 
two-quasiparticle strengths are shown in the lower panel. The arrow beside the 
abscissa axis indicates the neutron threshold energy 2|A| = 2.31 MeV. Right: 
Two-quasiparticle excitations generating the low-lying K'^ = 0+ modes at 3.3 MeV 
and 3.9 MeV. 

Table 3 



QRPA amplitudes of the K"" = 0+ mode at 3.3 MeV in ^SMg. This mode has 
B{E2) = 1.67 e^fni^ B(Q^2)=66.3 fm^, B{Q^^2)=89.0 fm^, and Yl l5a/3p = 2.32 x 
10~^. Only components with l/a/sp — |5a/3p > 0.01 are listed. 





a 


/3 


(MeV) 


If |2 1^ |2 
\ja/3\ - WafSl 


(fm2) 


(fm2) 


(a) 


z/[310]l/2 


i/[310]l/2 


3.37 


0.673 


6.08 


5.25 


(b) 


z/[312]5/2 


i/[312]5/2 


4.84 


0.146 


0.821 


-0.293 


(c) 


z/[310]l/2 


1/ [330] 1/2 


5.35 


0.023 


-3.59 


0.769 


(d) 


z/[303]7/2 


i/[303]7/2 


6.35 


0.066 


-2.64 


0.614 


(e) 


i/[202]3/2 


i/[202]3/2 


7.82 


0.021 


-1.29 


0.149 



portance of dynamical pairing effects in generating soft dipole excitations has 
been demonstrated by Matsuo et al. [13] for spherical unstable nuclei near the 
neutron drip line. 



Next we discuss the = O"*" excitations in "^^Mg and ^^Mg. The quadrupole 
transition strengths calculated for ^^Mg arc presented in Fig. 5, which exhibits 
two peaks below 4 MeV. The major two-quasiparticle excitations generating 
these peaks are illustrated in the middle and right panels of this figure. Their 
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Table 4 

QRPA amplitudes of the K"" = 0+ mode at 3.9 MeV in ^'^Mg. This mode has 
B{E2) = 4.72 e^fm^, B{Q''2)=68.1 fm'^, B{Q^^2)=m fm^ and Yl lfa/3p = 2.71 x 
10^^. Only compovunils with l/a/jp — l^a/jp > 0.01 are listed. 





a 


/3 


(MeV) 




(fm2) 


2U,ap 

(fm2) 


(a) 


i/[310]l/2 


i/[310]l/2 


3.37 


0.037 


6.08 


1.34 


(b) 


i/[301]l/2 


i/[310]l/2 


4.42 


0.258 


1.67 


-1.20 


(c) 


z/[312]3/2 


i/[312]3/2 


4.90 


0.048 


0.716 


0.169 


(d) 


z/[312]3/2 


i/[321]3/2 


5.47 


0.250 


-3.04 


-2.20 


(e) 


z/[301]l/2 


i/[301]l/2 


5.47 


0.018 


0.802 


0.131 


(f) 


z/[321]3/2 


i/[321]3/2 


6.04 


0.058 


1.66 


-0.411 


(g) 


z/[303]7/2 


z/[303]7/2 


6.35 


0.084 


-2.64 


-0.853 


(h) 


i/[330]l/2 


i/[330]l/2 


7.33 


0.099 


4.57 


-1.48 



Table 5 

QRPA amplitudes of the K'^ = 0+ mode at 3.4 MeV in ^^Mg. This mode has 
B{E2) = 8.1 e2fm^ B{Q''2)=104 fm^ B{Q^^2)=170 fm*, and Eba/3p = 3-91 x 
10~^. Only components with l/a/sP — Iffa/sP > 0.01 are listed. 





a 


P 


(MeV) 




(fm2) 


(fm2) 


(a) 


z/[310]l/2 


z/[310]l/2 


4.06 


0.071 


5.80 


-1.58 


(b) 


z/[321]3/2 


z/[321]3/2 


4.48 


0.098 


4.60 


-1.61 


(c) 


u[S12]5/2 


z/[312]5/2 


4.87 


0.227 


0.714 


0.347 


(d) 


i/[310]l/2 


1/ [330] 1/2 


4.91 


0.211 


-3.08 


-2.11 


(e) 


i/[301]l/2 


i/[310]l/2 


5.69 


0.033 


2.02 


-0.511 


(f) 


i/[330]l/2 


1/ [330] 1/2 


5.76 


0.116 


3.98 


-1.50 


(g) 


[202] 3/2 


[202] 3/2 


5.79 


0.046 


-1.47 


-0.271 


(h) 


[303] 7/2 


[303] 7/2 


7.67 


0.049 


-1.82 


-0.411 


(i) 


7r[211]l/2 


TT [220] 1/2 


6.44 


0.054 


-0.251 


-0.599 



QRPA amplitudes are listed in Tables 3 and 4. From these Tables, it is seen 

that the peak at 3.3 MeV is mainly generated by the particle-particle type 
z/[310]l/2®z/[310]l/2 and z/[312]5/2®z/[312]5/2 excitations, while the peak at 
3.9 MeV is mainly associated with the particle-hole type i/[301]l/2®i/[310]l/2 
and i/[312]3/2 ® z/[321]3/2 excitations. 
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The quadrupole transition strengths calculated for '^^Mg are displayed in 
Fig. 6. We notice a prominent peak at about 3.4 MeV below the one-neutron 
threshold energy (4.1 MeV), which possesses a strongly enhanced transition 
strength of about 24 W.u. (1 W.u. ~ 7.1 fm^ for ^^Mg). This peak exhibits 
a clear character of collective vibration: As seen from Table 5, this collective 
mode is created by coherent neutron excitations. Its main components arc the 
particle-hole type z/[310]l/2®z/[330]l/2 and z/[301]l/2®z/[310]l/2 excitations 
and the particle-particle type i/[312]5/2(8)i/[312]5/2 and i/[321]3/2(g)i/[321]3/2 
excitations. These particle-particle type and particle-hole type excitations are 
coherently superposed to generate this collective neutron mode. 



3.3 K"" = 2+ modes 



Let us now turn to the K'^ = excitation modes. The quadrupole transition 
strengths calculated for ^°Mg are displayed in Fig. 7. We notice a prominent 
peak at about 2.8 MeV which possesses strongly enhanced transition strength 
of about 19 W.u. The QRPA amphtudes of this excitation are listed in Table 6. 
Prom this Table, we see that this peak represents a collective excitation consist- 
ing of a coherent superposition of the proton particle-hole excitation from the 
[211]3/2 level to the [211]l/2 level and a number of neutron two-quasiparticle 



180 
150 
120 
_90 
E 60 
^ 30 

-I—' 

^ 

S50 
w 

40 
30 
20 
10 




isoscalar ^^^g "^^ 



-18 



unperturbed 



% -20 



c 



1 2 3 4 5 
hco (MeV) 



-22 



-24 



-26 



[211]1/2 



[220] 1/2 



-1 - 

-2 

-3 



[301] 1/2 
[303]7/2 ^ 



[310]1/2 
[312]5/2 



[321]3/2 



[330] 1/2 
[202]3/2 



proton 



neutron 



Fig. 6. Left: Isoscalar quadrupole transition strengths B{Q^^2) for the K'^ = 0"*" exci- 
tations in ^^Ivlg are plotted in the upper panel, while unperturbed two-quasiparticle 
strengths are shown in the lower panel. The arrows beside the abscissa axis indicate 
the neutron threshold energy £^th = 4.06 IVleV (Iqp continuum) and 4.12 MeV (2qp 
continuum). Right: Two-quasiparticle excitations generating the lowest K'^ = 0+ 
mode at 3.4 MeV in ^^Ivlg. 
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Fig. 7. Left: Isoscalar quadrupole transition strengths B{Q^^2) for the K'^ = 2+ 
excitations in ^^Mg. Results of the QRPA calculation with and without including 

the dynamical pairing effects are plotted in the upper and middle panels, respec- 
tively, while unperturbed two-quasiparticle strengths are shown in the lower panel. 
Notice that different scale is used for the unperturbed strengths. The arrow beside 
the abscissa axis indicates the neutron threshold energy 2|A| = 0.82 MeV. Right 
Two-quasiparticle excitations generating the lowest K'^ = 2^^ mode at 2.9 MeV. 
Two-quasiparticle excitations satisfying the asymptotic selection rule for the 7 vi- 
bration (AA^ = 0, An3 = 0, AA = 2) are drawn by thick arrows. 

excitations. Similarly to the K'^ = O"*" excitation modes discussed in the pre- 
vious subsection, the asymptotic selection rule (AA^ = 0, Ana = 0, AA = 2) 
well known for the 7 vibrations is violated for the neutron excitations, because 
these levels are loosely bound or resonances and their quasiparticle wave func- 
tions are significantly extended outside of the nucleus. On the other hand, pro- 
ton particle-hole excitations satisfy the selection rule because they are deeply 
bound. Wc also show in Fig. 7 the result of QRPA calculation where the resid- 
ual pairing interaction is turned off. Comparing with the full QRPA result, we 
see that the transition strength is reduced about 30%. Thus, the dynamical 
pairing effect is important, though its effect is weaker than for the K'^ — 0+ 
mode. This is because the K'^ = 2"*" mode consists of both proton and neutron 
excitations and the pairing is effective only for neutrons. 

The quadrupole transition strengths calculated for ^^Mg and ^^Mg are dis- 
played in Figs. 8 and 9, respectively. For each case, we see a prominent peak at 
about 2.9 MeV which possesses strongly enhanced transition strength (about 
15 W.u. and 12 W.u. for '^^Alg and "^^Mg, respectively). The QRPA amplitudes 
of these modes are listed in Tables 7 and 8. These modes possess essentially 
the same microscopic structure as the collective K'^ — 2+ mode in ^°Mg dis- 
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Table 6 

QRPA amplitudes of the K'^ = 2+ mode at 2.9 MeV in ^°Mg. This mode has 
B{E2) = 11.7 e^fm^, B{Q''2)=75.7 fm"^, B{Q^^2)=U7 fm"^, and Yl l5a/3p = 6.73 x 
10"^. Only components with |/a/3p - |5a/3p > 0.01 are listed. The label 2/1/2" 
denotes a discretized non-resonant continuum state. 





a 




(MeV) 


/a/3 P - 5a/3 P 


'^22,a,i3 

(fm2) 


22,ap 

(fm2) 


(a) 


;/[312]3/2 


;/[310]l/2 


3.77 


0.013 


1.22 


-0.145 


(b) 


z/[301]l/2 


i/[312]3/2 


4.16 


0.098 


-5.37 


-1.75 


(c) 


z/[310]l/2 


i/[312]5/2 


4.51 


0.085 


-4.37 


-1.34 


(d) 


z/[312]3/2 


i/[303]7/2 


4.88 


0.011 


-5.03 


-0.454 


(e) 


z/[301]l/2 


i/[312]5/2 


4.90 


0.016 


-2.07 


-0.296 


(f) 


i/[310]l/2 


i/[321]3/2 


5.34 


0.047 


-2.67 


-0.663 


(g) 


1/1/2- 


i/[312]5/2 


6.96 


0.015 


1.93 


-0.298 


(h) 


1/1/2" 


i/[321]3/2 


7.28 


0.018 


1.46 


-0.265 


(i) 


7r[211]l/2 


7r[211]3/2 


4.32 


0.596 


-2.11 


-2.02 
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Fig. 8. Left: Isoscalar quadrupole transition strengths B{Q^^2) for the K"^ = 2^ exci- 
tations in ^^Mg are plotted in the upper panel, while unperturbed two-quasiparticle 
strengths are shown in the lower panel. The arrow beside the abscissa axis indicates 
the neutron threshold energy 2|A| = 2.31 MeV. Right: Two-quasiparticle excitations 
generating the lowest K'^ = 2+ mode at 2.9 MeV in ^^Mg. 
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Table 7 

QRPA amplitudes of the K"" = 2+ mode at 2.9 MeV in ^'^Mg. This mode has 
B{E2) = 11.2 e^fm^, B{Q''2)=53.0 fm^, B{Q^^2)=n3 fm^ and ^ |c/a/3p = 6.95 x 
10^^. Only components with |/a/3p — ba/3p > 0.01 are listed. 





a 


P 


(MeV) 




(fm2) 


22,ap 

(fm2) 


(a) 


i/[310]l/2 


z^[312]5/2 


4.10 


0.146 


-3.89 


-1.54 


(b) 


i/[312]3/2 


jv[310]l/2 


4.13 


0.016 


-0.221 


-0.032 


(c) 


z/[310]l/2 


i/[321]3/2 


4.70 


0.108 


2.81 


-1.05 


(d) 


z/[301]l/2 


i/[312]5/2 


5.15 


0.010 


1.81 


-0.204 


(e) 


z/[301]l/2 


i/[312]3/2 


5.18 


0.031 


-3.58 


-0.690 


(f) 


z/[312]5/2 


1/ [330] 1/2 


6.08 


0.017 


-1.42 


-0.239 


(g) 


i/[321]3/2 


[330] 1/2 


6.68 


0.029 


1.61 


-0.344 


(h) 


7r[211]l/2 


7r[211]3/2 


4.43 


0.541 


-2.09 


-1.94 



Table 8 

QRPA amplitudes of the K"" = 2+ mode at 2.9 MeV in ^"^Mg. This mode has 
B{E2) = 10.6 e^fm^, B(Q''2)=35.3 fm^, B{Q^^2)=8'i.6 fm^, and ^ l^^gp = 7.06 x 
10~^. Only components with j/a/jp — baflP > 0.01 arc listed. The columns (f) 
and (f) are assigned the same configuration, because we obtain two discretized 
continuum states associated with the [200] 1/2 level for which > |A . 





a 




(McV) 


/q/3 P - ba/? P 


(fm2) 


^^^22,al3 

(fm2) 


(a) 


i/[310]l/2 


i/[321]3/2 


4.27 


0.087 


1.68 


0.552 


(b) 


z/[310]l/2 


z/[312]5/2 


4.46 


0.045 


-2.62 


0.551 


(c) 


z/[321]3/2 


1/ [330] 1/2 


5.12 


0.165 


2.51 


1.20 


(d) 


z/[312]5/2 


1/ [330] 1/2 


5.31 


0.034 


-1.36 


0.305 


(e) 


z/[301]l/2 


z/[312]5/2 


6.09 


0.011 


1.66 


0.202 


(f) 


i/[202]3/2 


1/ [200] 1/2 


6.90 


0.010 


-2.633 


0.270 


(f) 


i/[202]3/2 


1/ [200] 1/2 


7.23 


0.013 


0.981 


0.136 


(g) 


7r[211]l/2 


7r[211]3/2 


4.54 


0.547 


-2.09 


1.87 



cussed above. They also correspond to the 7 vibrational mode obtained in 
the previous QRPA calculation [32] for ^^Mg. In our calculation, however, the 
collectivity of these modes remains almost the same even if we use different 
deformations for protons and neutrons, differently from Ref. [32]. 
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Fig. 9. Left: Isoscalar quadrupole transition strengths B{Q^^2) for the 
K'^ = 2+ excitations in ^^Mg are plotted in the upper panel, while unperturbed 
two-quasiparticle strengths are shown in the lower panel. The arrows beside the 
abscissa axis indicate the threshold energies 4.06 McV (Iqp continuum) and 4.12 
MeV (2qp continuum). Right: Two-quasiparticle excitations generating the lowest 
K'^ = 2+ mode at 2.9 MeV. 



3.4 Dependence on pairing interaction 



In this subsection, we examine sensitivity of the low-frequency = 0+ and 2+ 
modes on the density dependence of the pairing interaction. For this purpose, 
we repeated the HFB and QRPA calculations using pairing interactions with 
density dependence different from the surface type (r^ = 1.0 in Eq. (2)); i.e., 
the mixed type {ri = 0.5) and the volume type (rj = 0.0). Since the result for 
the mixed-type pairing is intermediate between those for the surface-type and 
the volume-type, we show in Fig. 10 only the quadrupole transition strengths 
obtained using the volume-type pairing interaction. In this calculation, the 
pairing interaction strength Vq — —215.0 MeV-fm^ is chosen to yield approxi- 
mately the same average pairing gaps as those for the surface type. Comparing 
with the results obtained using the surface-type pairing, shown in Fig. 3, we 
see that the transition strengths for the K'^ = 0^ collective modes are appre- 
ciably reduced, while those for the K'^ = 2"*" collective modes are almost the 
same. We have checked that, although the strengths are reduced, the micro- 
scopic structure of these collective modes are basically the same as discussed 
above on the basis of the results of calculation using the surface-type pairing 
interaction. Thus, we can say that the quadrupole transition strengths for the 
low-frequency K'" = 0"*" collective modes are especially sensitive to the density 
dependence of the pairing interaction. Such a sensitivity has been stressed 
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Fig. 10. Isoscalar quadrupole transition strengths B{Q^^2) for the K = 0^ excita- 
tions (upper panel) and the if = 2+ excitations (lower panel) built on the prolately 
deformed ground states of 36,38,40]y[g r^^ie QRPA calculations are made in the same 
way as in Fig. 2, except that the volume-type pairing interaction is used here. The 
arrows indicate the neutron threshold energies; 3.69 MeV (Iqp continuum) and 3.71 
MeV (2cip continuum) for ^^Mg, 1.77 MeV (2qp continuum) for ^^Mg, and 0.15 MeV 
(2qp continuum) for ^°Mg. 

also by Matsuo et al. in their continuum QRPA calculations for El strength 
functions in neutron rich O, Ca and Ni isotopes[13]. 



4 Concluding remcirks 



We have carried out the QRPA calculations on the basis of the deformed 
WS plus HFB mean field in the coordinate representation, and obtained the 
low-frequency K"^ = 0+ and 2+ collective modes in deformed 36,38,40]y[g gjQgg 
to the neutron drip line. It has been shown that these modes possess very 
strong isoscalar quadrupole transition strengths. One of the reasons of this 
enhancement is that the quasiparticle wave functions participating in these 
collective excitations have spatially extended structure. The other reason is 
that the residual pairing interactions, in addition to the particle-hole type 
residual interactions, enhance the collectivity of these modes. The result of 
the present calculation suggests that the low-frequency — 0"*" collective 
mode is a particularly sensitive indicator of the nature of pairing correlations 
in nuclei close to the neutron drip line. 

This paper should be regarded as an exploratory work toward understand- 
ing low-frequency collective modes of excitation in unstable nuclei close to 
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the neutron drip line. It is certainly desirable to improve the treatment of 
the contimmm at least in the following points. First, one may try to use a 
smaller mesh size and a larger box by implementing an adaptive coordinate 
method[24]. Second, one may try to take into account the width of resonance 
by employing Gamow states as basis of the QRPA calculation [56]. The result 
of the present work indicates that calculations using such an improved frame- 
work will be very interesting and worthwhile. We plan to attack this subject 
in future. 
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A Eigenphase sum for single-particle resonance states 

We examine properties of three single-particle states in the continuum, which 
play a key role in generating the low-lying excitations modes in 36,38,40]y[g^ 
The resonance energy and width in a deformed potential can be estimated 
using the eigenphase sum A{E). It is defined in terms of the eigenvalues of 
the scattering matrix (S-matrix) as 

a 

We evaluate the eigenphase sum for three states following the procedure of 
Ref. [52]. The resonance energy and width are identified with the peak energy 
of ^dA{E)/dE and its FWHM, respectively [53,54]. This evaluation IS m good 
correspondence with another definition of the resonance; the Gamow state in 
a deformed potential [55] which represents the pole of the (S— matrix in the 
complex momentum plane. 

The result of this calculation, presented in Fig.A.l, indicates that the [301] 1/2 
and [312] 3/2 states can be regarded as resonances with rather large widths; 
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Fig. A.l. The eigenphase sum (upper panel) and its derivative (lower panel) for 
the il'^ = l/2~,3/2~ and 7/2~ states in "^°Mg are plotted as functions of energy. 

their energies are 0.53 — iO.46 (MeV) and 0.42 — iO.33 (MeV), respectively. 
On the other hand, the [303] 7/2 state is evaluated as a narrow resonance with 
energy 0.44 — iO.0005 (MeV). Obviously, the small width is due to its high 
centrifugal barrier. 
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